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A numerical renormalization-group survey of the zero-bias electrical conductance through a quantum dot 
embedded in the conduction path of a nanodevice is reported. The results are examined in the light of a recently 
derived linear mapping between the temperature-dependent conductance and the universal function describing 
the conductance for the symmetric Anderson model. A gate potential applied to the conduction electrons is 
known to change markedly the transport properties of a quantum dot side-coupled to the conduction path; in 
the embedded geometry here discussed, a similar potential is shown to affect only quantitatively the tempera- 
ture dependence of the conductance. As expected, in the Kondo regime the numerical results are in excellent 
agreement with the mapped conductances. In the mixed-valence regime, the mapping describes accurately the 
low-temperature tail of the conductance. The mapping is shown to provide a unified view of conduction in the 
single-electron transistor. 



I. INTRODUCTION 

The development of the first single-electron transis- 
tor (SET) 1 was preceded by analytical and numerical 
breakthrough s - 2 ' 3 ' 4 ' 5,6 ' 7 and, in the subsequent years, motivated 
numerous theoretical investigations ! 8 ' 9 ' 10 ' 1 V 2 ' 13,14 ' 15 ' 16 ' 17 ' 18 
Notwithstanding this intense activity, important aspects 
of the transport properties of the device received limited 
attention. In particular, even after it guided the interpretation 
of experimental data drawn out of nanodeviceS j 19,20 - 21 ' 22 the 
concept of universality remained confined to the narrow 
corner in which it was first identifiedi 2 - 

As it was established much more recently, in the Kondo 
regime of the Anderson model, 23 the thermal dependence of 
the zero-bias conductance of quantum dots either embedded 24 
in or side-coupled 25 to the conduction path of elemental nan- 
odevices maps linearly onto the conductance curve for the 
symmetric model computed by Costi, Hewson and Zlatic. 2 
The linear coefficient in the mapping, which is parametrized 
by the ground-state phase shift 6 of the conduction electrons to 
which the quantum dot is coupled, depends on the geometry. 
Qualitatively different thermal dependences result, which dis- 
tinguish the side-coupled geometry from the embedded con- 
figuration. With a side-coupled quantum dot, the linear co- 
efficient depends only on the phase shift. As it results, the 
application of a potential W to the conduction electrons can 
switch the conductance curves from monotonically increasing 
to monotonically decreasing functions of the temperature. 25 

In the embedded configuration, of which the SET is the 
simplest illustration, the conductance always decreases with 
temperature. The mapping depends on the difference 6 - 6 W 
between the ground-state phase shift and the phase shift 5 W 
that the potential W would induce if the conduction electrons 
were decoupled from the dot. In the special class of model 
Hamiltonians that are invariant under particle-hole transfor- 
mations, the symmetry makes the dot occupation tij unitary, 



and it follows from the Friedel sum rule* 6 - that 6 — n/2. For 
asymmetric Hamiltonians, by contrast, a priori knowledge of 
6 is impossible. In the Kondo regime, again on the basis 
of the Friedel sum rule, a difference 6 - 6 W close to n/2 is 
expected; 24 accurate estimates nevertheless require diagonal- 
ization of the model Hamiltonian. One possibility would be 
to generate temperature-dependent conductance curves from 
Bethe-ansat z 27 ' 28 results for the phase shift. Here, however, 
we prefer the numerical renormalization-group (NRG) ap- 
proach, which gives direct access to the phase shifts and to 
the thermal dependence of the conductance. 

This paper compares numerically computed SET conduc- 
tances with the mapping to the universal function. From 
the same NRG diagonalization of the model Hamiltonian, we 
compute (i) the conductance as a function of the temperature; 
and (ii) the ground-state phase shift, which specifies the map- 
ping. Plotted as functions of the temperature T in the Kondo 
regime, the computed conductances rise from from nearly 
zero to nearly ballistic as T decreases past the Kondo temper- 
ature, and conduction electrons screen the the dot magnetic 
moment. The curves obtained from the mapping run through 
the numerical data. 

Outside the Kondo domain, the mapping to the univeral 
function is expected to fail. To witness its downfall, the nu- 
merical survey includes the adjacent mixed- valence domain, 
in which the dot moment is only partially formed. Here, the 
conductance crosses over to its ground-state value at relatively 
high temperatures, which are outside the domain of the map- 
ping; hence, only the low-temperature tail of the numerical 
data can be accurately fitted. 

Our presentation is distributed over five Sections. SectionHIl 
describes the SET and the Anderson Hamiltonian modeling it. 
Section|III]discusses cursorily the mapping derived in Ref.l24l 
Section[lV]is dedicated to the NRG procedure, which it sum- 
marizes, and to the numerics, which it details. The numerical 
results appear next, in Section [V] which starts with a discus- 
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sion of the phase-shift differences 6 - 6 K . 

The same section presents the conductance curves. Be- 
sides describing quantitatively the conductance in the Kondo 
regime (and the low-temperature sector of the mixed-valence 
regime), the mapping to the universal function offers a simple, 
unifying view of charge transport through a single-electron 
transistor. Faithful to this notion, Sections fV Al and TV Bl com- 
pound the output of 100 NRG runs in two plots of the con- 
ductance as a function of the temperature and gate potential 
applied to the quantum dot; and show that every feature of the 
landscapes is easily understood in the light of that mapping. 
Section lVCl then turns the numerically computed ground-state 
phase shifts into linear coefficients, and the mapping into con- 
ductance curves that fit the NRG results for the temperature- 
dependent conductances. Our conclusions and a summary 
constitute the closing Section IVTI 



H. MODEL 

Figure Q] depicts the object of our study, a quantum dot 
symmetrically coupled to two otherwise independent quantum 
wires. The tunneling amplitude V transfer charge between the 
dot and the wires. To represent the dot, we introduce a sin- 
gle, spin-degenerate level c ,. An energy s d , controlled by the 
gate potential Vd, and a Coulomb repulsion U define the dot 
Hamiltonian Hd, which can be written 35 



U 



U 



H d = (Sd + -j)nd - -jiridi - n d[ ) , 



(1) 



to emphasize that the energy s d + U/2 breaks the particle-hole 
symmetry of the dot Hamiltonian. 
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FIG. 1 : (Color online) Single-electron transistor. The gate potentials 
V g and V v control the dot energy sj and the energy of the Wannier or- 
bital /„ , respectively. Symmetric tunneling amplitudes V couple the 
dot to the quantum wires and allow conduction along the indicated 
path. 



and W is the potential due to the wire electrodes in Fig.Q] 

The coupling V, to the wires, broadens the dot level c d . The 
scattering potential W reduces its width F = npV 2 , where Np 
is the density of conduction states, t< 
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r„ = 



1 +7T 2 p 2 W 2 ' 



(4) 



With W — 0, a specially important instance of Eq. (f2| is the 
particle-hole symmetric Hamiltonian 

Hl=Y A ^ala k +Wflf +V(jlc d +n. C ^-j^-n l f, (5) 



to which Ha reduces for s d + U/2 = W — 0. 

The operators a k on the right-hand side of Eq. ((2j are even 
combinations of the conduction operators in the left and in the 
right wire segments in Fig.[T] a k — (ql + qr)/ y2. The inver- 
sion symmetry of the device decouples the odd combinations 
(ckL—CkR)/ V2 from the dot level, which make no explicit con- 
tribution to the zero-bias conductance G(T). Linear Response 
relates G(T) to the dot-level spectral density p d ? A 



G(T) = Q 2 



nV K f 

J-D 



Pd(e, T) 



dm 

de 



de, 



(6) 



where Q2 = 2e 2 /h is the quantum conductance, /(e) is the 
Fermi function, Z is the partition function for the Hamiltonian 
Ha, and 



Pd(e, T) = 1 J] e -jrxM c\ \m)\ 2 5{E m - E„ - e). (7) 



Here, \m) (\n)) denotes an eigenstate of Ha with eigenvalue 

Em (En). 

The substitution of Eq. (jT) on the right-hand side of Eq. (0 
yields an expression that translates into a few lines of com- 
puter code: 



r(T , r P*r w v \{m\c d \n)\ 2 
G(T) =§2— 2, ^7 



z 



e fSE„ ■ 



(8) 



As this expression suggests, the computational efforts under- 
lying a conductance curve G(T) and e. g., a magnetic suscep- 
tibility plot 29 are comparable. 



Inclusion of the quantum wires, represented by a structure- 
less half-filled band of width 2D containing conduction 
states, and of their coupling to the quantum dot yields the An- 
derson Hamiltonian 

H A = J] e k a\a k + W$f Q + V($c d + H. c.) + H d , (2) 

k 

where 

/ =£a*/ViV, (3) 

k 



III. MAPPING 

Equation ([8]) yields conductances ranging from zero to the 
quantum conductance We are interested in the Kondo 
regime, the vast parametric subspace defined by the condi- 
tion F w <K minde^l, U + Sd) and ksT «: min(|e^|, U + Sd,D). 
At the high-temperature end of the Kondo regime a magnetic 
moment arises, associated with the nearly unitary dot-level oc- 
cupancy Hd- As the temperature is lowered past the Kondo 
temperature T K , and the wire electrons screen the resulting 
dot magnetic moment, the Kondo cloud couples the dot level 
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to the wire states and sustains conductance. Costi, Hewson 
and Zlatic 2 showed that the thermal dependence of the con- 
ductance for the symmetric Hamiltonian (0 is a universal 
function G S (T /T£) of the temperature scaled by the Kondo 
temperature. More recently, we have shown that the Kondo- 
regime conductance maps linearly onto G s (T /Tk)'M 

c (£)-f=-( cS £)-fH"-«- <" 

where S is the ground-state conduction-band phase shift, and 
6 W , the Fermi-level phase shift for V = 0, i. e., 



tanc) M , = -npW. 



(10) 



As the temperature rises, the universal function decays from 
G S (T « T K ) = g 2 , through G S (T = T K ) = g 2 /2, to 
G S (T » Tk) = 0. The linear coefficient on the right-hand 
side of Eq. © is a function of the phase-shift difference 6-6 w . 
It follows from the Friedel sum rule that 2(6 - 6 W ) - n d 
so that in the Kondo regime, and even in the mixed-valence 
regime, S - S w is never small. According to Eq. (0, for 
6 - 5 W = tt/2, the conductance G(T/Tk) sticks to G s (T/Tk)- 
A particular case is the symmetric Hamiltonian (0, for which 
8 = tt/2, 6 w = 0, and G(T/T K ) = G S (T/T K ). If the differ- 
ence were 5 - 6 W - n/4, on the other hand, the conductance 
would be flat: G(T) = Q 2 I2. For the intermediate differences 
7r/4 < d- 5 W < n/2 observed in the Kondo and mixed-valence 
regimes, the conductance is a monotonically decreasing func- 
tion of the temperature. If 5-5 w + n/2, G(T <sz Tk) is smaller 
than the quantum conductance, G(T » 7V) is nonzero, and 
G(T/T K ) is flatter than G S (T/T K ). 

Equation (O offers a qualitative view of the thermal depen- 
dence of the conductance, which combined with the Friedel 
sum rule, describes G(T/Tk) semiquantitativelyi^i A more at- 
tractive alternative, the numerical evaluation of 6 and Tk is 
discussed next. 



IV. NUMERICS 

Excellent descriptions of the NRG method being 
available ^ 8 ' 29 ' 30 brief recapitulation of the four steps consti- 
tuting the procedure will be sufficient. Two dimensionless 
parameters A > 1 and < z < 1 define the logarithmic 
discretization of the conduction band£2J. The infinite energy 
sequence E,„ = DA l ~ z ~ m (to = 0,1,...) defines the inter- 
vals I m = [E m+ i,E,„]. For each interval, a single operator 
a m + - pjjj cikdekln m , with normalization factor n m , is 
defined. In the negative half of the conduction band, the 
sequence — E m (to — 0, 1 . . .) defines the mirror image a,„_ of 
each operator a m+ . The a m± forma a basis upon which the 
conduction band Hamiltonian is projected! 2 ^ 

Next, a Lanczos transformation 32 makes tridiagonal the 
projected conduction Hamiltonian, so that the model Hamil- 
tonian reads 

Ha = If] tnfnUi + yf c d + H. c. + Wflf + H d . (11) 



Here, / is the operator defined in Eq. (O, and the f n 's 
(n = 0, 1 , . . .) form an orthonormal basis that replaces the 
a m± 's (m = 0, 1, . . .). With z = 1, we recover the Lanczos 
transformation in Ref. |25j. Otherwise, the codiagonal coef- 
ficients t n have to be determined numerically. 31 With error 
0(A-"), it is founds that t n = D[(l - A" 1 )/ log A] A 1 ^-" 12 . 
This shows that the f„'s decrease rapidly with «, a conclusion 
that brings us to the third step in the NRG procedure, the def- 
inition of a renormalization-group transformation. 

Given a temperature T and a small dimensionless parame- 
ter a, let be the smallest integer such that t N < akgT, and 
consider the infinite sum on the right-hand side of Eq. ( fTTb . 
Compared to k^T ', the codiagonal element t^ is then negli- 
gible, and to compute G(T), it is safe to neglect the term 
with n - N. This decouples the subsequent terms in the 
sum from the quantum dot, so that they no longer contribute 
to the conductance. The infinite sum can therefore be trun- 
cated at n = N - 1. We define the reduced bandwidth 
D N = D[(l - A _1 )/logA]A- (JV - 1)/2 and the dimensionless, 
scaled, truncated Hamiltonian H% : 



D " H % = Z tnf » fn+1 + Vf o c d + H - c - 



-Wflf )+ H d . (12) 



The entire Hamiltonian has been divided by so that the 
smallest codiagonal coefficient in the scaled sum, t^-\ID^, is 
of 0(1). 

The last step in the NRG procedure is the iterative diago- 
nalization of the model Hamiltonian. With N = 0, the right- 
hand side of Eq. (fT2l) is easily diagonalized; four eigenval- 
ues and four eigenvectors |m) (to = 1, ... ,4) result. At 
this stage, it is equally simple to calculate the matrix elements 
(m] c d |n) between the eigenvectors of H^ =0 , which will be 
needed to compute the right-hand side of Eq. ([8]). 

Application of the operators f^, fl, f^f^, and the iden- 
tity 11 on the eigenvectors of H? =0 generates sixteen states that 
constitute a basis upon which the Hamiltonian H^ =l can be 
projected. Appropriately chosen linear combinations of those 
operators yield basis states \p)\ (p — 1, . . . , 16) that diagonal- 
ize the charge and spin operators; projected on them, H^ =1 
reduces to block-diagonal matrices, which are then diagonal- 
ized numerically. The matrix elements Q (m\ c d \n) Q (to, n 

I - 1 ) are projected onto the basis \p)\ and subsequently 

rotated to the basis of the eigenstates \m)\ (to = 1, . . . , 16) 
of H^ =l . Application of the operators fL, f^f^, and 11 
on the |to)i creates 64 basis vectors upon which H^ =2 can be 
projected, and the procedure is iterated. 

To check the exponential growth of matrix dimensions, a 
dimensionless parameter v is chosen, which will control the 
cost and the accuracy of the iterative diagonalization. At 
the end of iteration N, the eigenvectors with scaled energies 
Em/DN above v are discarded before the construction of the 
basis states \p)n+i, upon which the Hamiltonian H^ +l will 
be projected. This expedient limits the number of basis states 
and hence the computational effort that each iteration requires. 
The cost of a full NRG run grows linearly with the number of 
iterations. 
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The diagonalization yields scaled eigenvalues ranging from 
unity to v, i. e., a window of energies ranging from D N to vD N . 
Having neglected t N m D^+i, we can only compute conduc- 
tances for k B T > yDtf + \, where y > 10. At the other extreme, 
the ultraviolet truncation restricts us to temperatures such that 
k B T < vD N . Thus, provided that (D N v)/(yD N+ i) > VX, 
i. e., that v > y, the A-th iteration yields reliable conduc- 
tances in the temperature window yD^+i < k B T < ^fKyD^ + \. 
If a run is stopped at iteration N max , the juxtaposition of the 
resulting windows yields G(T) for all temperatures above 
yD Nmax+ i/k B . In practice, the conductance is only computed 
for the k B T < 0.1 D, because irrelevant operators artificially 
introduced by the logarithmic discretization make the interval 
0.1 D < k B T < D unreliable. 



Conductance curves computed with large A show oscil- 
lations, which can be traced to a sequence of poles on the 
3e = ±;7r/logA lines of the complex-energy planed To 
eliminate these artifacts of the discretization, we average the 
conductance curve G(T) computed for given z over a sequence 
of equally spaced z's in the interval < z < l. 34 The expo- 
nential deependence of the computational effort on 1 / log A, 
makes this averaging procedure far more efficient than com- 
parably accurate computations with small A. 



The conductances in SectionlVlwere computed with A = 6 
and averaged over two z's: 0.5 and 1. The amplitude of the 
residual oscillations encountered after averaging over z, some- 
what smaller than 0.001 e 2 /h, provides an estimate of the er- 
ror introduced by the logarithmic discretization. The other 
two parameters controlling the precision of the results were 
fixed at y = 10.5 and v = 50, respectively. Spin degenera- 
cies not counted, the number of states below the cutoff in each 
iteration peaked at 4000 in iteration 6. To estimate the er- 
ror due to the infrared and ultraviolet truncations, for each 
N < Nmax we compared the conductance at the lowest tem- 
perature in the (N - l)-th window, yD^ < k B T < VXyDw, 
with the conductance at the highest temperature in the A-th 
window. The mismatch between the two results never exceed- 
ing 0.001 e 2 /h, we conclude that deviations due to the three 
approximations in the procedure, the logarithmic discretiza- 
tion and the infrared and ultraviolet truncations, are compara- 
ble. At any temperature, the estimated absolute deviation in 
the computed conductances is smaller than than 0.05% of the 
quantum conductance. 



The relatively large discretization parameter expedites the 
calculation. On a standard desktop computer, a complete run, 
including (i) the iterative diagonalization of Ha and computa- 
tion of the matrix elements on the right-hand side of Eq. ^ 
for for each z, and (ii) the evaluation of the conductance curve 
in the interval 10~ 10 D < k B T < 0.1 D, takes less than 30 sec- 
onds. 



Phase shifts 

Along with the iterative diagonalization procedure, Eq. (fT2] i 
defines a renormalization-group transformation^ 

N+l 

T[<] = < +2 =AH N A+ Y-£- (/t/ n+1 + H. c.) . (13) 

The factor A multiplying the first term on the right-hand side 
magnifies the scale on which the eigenvalues of H% are ex- 
amined. On the new scale, the second term is a fine structure. 
In the absence of characteristic energies, as N grows the mag- 
nification compensates the refinement, and the lowest-energy 
eigenvalues of H^ +2 rapidly become indistinguishable from 
those of H%. This indicates that the Hamiltonian has reached 
a fixed point of T. 

In the Kondo regime, the condition V — turns the Ander- 
son Hamiltonian Ha into the local-moment fixed point (LM) 
of T. With V + 0, as the temperature is reduced past the domi- 
nant characteristic energy E c - min(|ej|, sj+U, D), the Hamil- 
tonian H^ first approaches the LM and then moves away to- 
wards the frozen-level fixed point (FL) — a strong-coupling 
fixed point equivalent to Eq. (f2]i with V — > oo. Between the 
LM and the FL lies the Kondo temperature T%, around which 
the conduction electrons screen the dot moment. 

If one of the dot excitation energies, Ao = or 
A2 = Ed + U, is smaller than the dot width T,,,, the model 
Hamiltonian enters the mixed-valenc regime. 38 Instead of 
min(Z), \sd\, \ed+U), the dominant characteristic energy is now 
E c = min(D, r vv ). The dot moment is only partially formed, 
as the coupling F w drives the model Hamiltonian toward the 
FL before it can come close to the LM. 

Devoid of characteristic energies, the two fixed points, LM 
and FL, are phase-shifted conduction bands, which can be di- 
agonalized analytically. 35 With V — 0, for instance, the model 
Hamiltonian flows towards the LM with phase shift 6 W . With 
V + 0, (inaccurate) estimates for the LM phase shifts can be 
extracted from the eigenvalues of Wl, where N is such that 
Dn » k B Tx, i. e., such that the model Hamiltonian dwells in 
the vicinity of the LM. 

Much more accurate FL phase shifts 6 can be obtained from 
the low-energy eigenvalues of H%, because for large N the 
eigenvalues of H% come arbitrarily close to the many-body 
excitations of the FL Hamiltonian 24 ' 35 

flk=5X*U*- (14) 

t,± 

Here, the + and - subscripts distinguish the positive eigen- 
values from the negative ones, while t = 0,1,... counts the 
positive (negative) eigenvalues upward (downward) from the 
Fermi level. 

Once the eigenvalues of H% are identified with the many- 
body energies generated from Eq. ( TPfl i. the ground-state phase 
shift 6 are extracted from the approximate expression describ- 
ing all but the e* (± closest to zero^ 5 - For A > 5, in particular, 
within 0.1% deviation, 

el = ±K V+Mln ({=1,2,.. .), (15) 



where v = 1— z(v=3/2~z) for odd (even) 2V. 



V. RESULTS 

To emulate the conditions under which a SET operates, we 
fix the Coulomb repulsion (U — 5 D) and effective dot-level 
width (T K - 0.15 D) and examine the ground-state phase shift 
and the thermal dependence of the conductance as a function 
of the dot energy Sd for five wire potentials W. Since Ha, 
G, and \6\ are invariant under the transformation cj — > — c' d , 

a k -> a\, s d + U/2 -> -(e d + U/2), W -> -W, we need 
not study negative wire potentials, which would mirror the 
conductances and phase shifts calculated with positive W. 

We turn first to the calculated phase shifts. Figure [2] dis- 
plays the argument of the trigonometric function on the right- 
hand side of the mapping (0, computed for five wire po- 
tentials W, in the dot-energy range < sj/D < -U. For 
W = 0, the upright triangles draw a well-known curve, 27 
which remains invariant under the particle-hole transforma- 
tion Sd + U/2 — > -(e,/ + U/2), 6 — > n — 6. At the symmet- 
ric point Sd + U/2 - 0, which corresponds to Eq. Q, the 
phase-shift is exactly n/2. The arrows above the top axis indi- 
cate the Kondo domain, within which the phase shift remains 
close to n/2. As \s d + U/2\ grows, the model Hamiltonian 
first approaches the limits of the Kondo domain and then in- 
vades the mixed-valence domain. In response, 5 moves away 
from n/2, towards zero for Sd + U/2 — > U /2, or towards n for 
s d + U/2 -> -U/2. 

The wire potential reduces the ground-state phase shift 
throughout the depicted range. For pW - 1, for instance, the 
phase shift at the symmetric dot-level energy Sd = -2.5 D — 
-U/2 is reduced from S - n/2 to 6 — n/10. In the Kondo- 
regime, as the illustration shows, the difference 6 - 6 H . is 
nonetheless pinned at n/2. The pinning is due to the Friedel 
sum rule. 26 Since the ground-state phase shift would be 6 W if 
V were zero, 2(6 - 6 w )/n is the screening charge due to the 
coupling to the dot. In the Kondo regime, that charge is nearly 
unitary, and 6 - 8 W « n/2. 

Figure [2] shows that a positive wire potential tends to dis- 
place the Kondo regime towards higher dot energies. For 
W — 0, the rapid decay of the phase shift near Sd = 
(sd = —U) marks the resonance between the rid — and 
rid — 1 (n d = 1 and rid = 2) dot-level configurations. The 
Kondo regime lies between them. As W grows, the two reso- 
nances move to higher e/s, and so does the Kondo regime. 



A. Conductance landscape 

According to Eq. Q, 5 - S w controls G(T). Consequently, 
the central features of Fig. [2] are manifest in landscape plots 
of the conductance. Figure [3] shows G(T) in the dot-energy 
range \e d + U/2\ < U/2 for U = 5 D, and T = 0.15 D. The 
plot surveys the entire Kondo regime and part of the mixed- 
valence regime. The plane e d - -U/2, which represents the 
symmetric Hamiltonian (0, splits the landspace in two sym- 



Kondo 



(W = 0) 

ed+U/2 




FIG. 2: (Color online) Ground-state phase shift S, measured from the 
phase shift <5„ obtained from Eq. JlOt for the displayed wire poten- 
tials W, as a function of the dot-level energy s d . The Friedel sum 
rule shows that the ordinate is equal to the dot occupation n c i- The <5's 
were obtained, with the help of Eq. d 1 5b , from the low-energy spec- 
trum of H A resulting from NRG runs with U = 5 D and T„. = 0. 15 O. 
The arrows above the top horizontal axis define the Kondo domain 
for W = 0. For W + 0, the Kondo domain is displaced to the right. 
Each vertical arrow pointing to the lower horizontal axis identifies 
the figure displaying the thermal dependence of the conductance for 
the indicated dot energy. 



metric halves, mapped onto each other by the particle-hole 
transformation Crf — > -ct, —> a\. 

Examined at the symmetric point Ed = -U/2, the 
temperature-dependent conductance reproduces the universal 
function G s (T/Tk). Here and elsewhere in the Kondo regime, 
the conductance at fixed Sd rises from zero to ballistic or 
nearly ballistic as the temperature is reduced past Tk, i- e., 
as one climbs from the high-temperature Coulomb-blockade 
valley to the low-temperature Kondo plateau. The Kondo tem- 
perature for the symmetric Hamiltonian is 8 x 10~ 7 D. As 
\sd + U/2\ grows, it rises until kgTk w 0.1 D, an equality indi- 
cating proximity to the mixed-valence regime, i. e., to the two 
resonances centereced at Sd = -5 D and s d = 0. As \sd + U/2\ 
grows further, we come into mixed-valence domain. The dot 
moment shrinks, and so does the Kondo cloud. The Kondo 
bypass of the Coulomb blockade becomes less and less effec- 
tive, and the conductance approaches zero. The steep drops 
near the s d = and s d = -5 planes in Fig.[3]mark the mixed- 
valence regime. 



B. Wire potential 

Figure [4] displays the conductance as a function of e d and 
T for U = 5D, F H , = 0.15 D, and pW = 1. Quantitative 
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FIG. 3: (Color online) Conductance as a function of the temperature 
and dot-level energy, for U = 5 D, T = 0.15 A and W = 0. The 
plot is symmetric with respect to the ej = — U/2 = -2.5 D plane. 
The sharp drops near e c i = -5 D and Sd = mark the borders of the 
Kondo regime, which extends roughly from sj = -T to sj + U = 
-T. In the Kondo regime, at fixed s d , the more gradual decay of 
the conductance with the temperature portrays the evaporation of the 
Kondo droplet. 
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FIG. 4: (Color online) Conductance as a function of the temperature 
and dot-level energy for U = 5 D, F„ = 0.15 D, and pW = 0.50. The 
wire potential breaks the particle-hole symmetry visible in Fig. [3] 
The sharp drop centered at ej = -5D in Fig.[3]is now fully visible, 
while the one centered at Ed = is out of sight, an indication that the 
Kondo regime has been displaced to higher dot energies. The bell- 
shaped resonance near the bottom left corner of the k B T = 10 1 D 
plane stakes the mixed-valence regime. 



differences distinguish the plot from Fig. [3] In particular, the 
Kondo temperature is now minimized at the higher dot-level 
energy sj = -1.9 D, the minimum kgTx - 2.4 x 10~ 6 D is 
thirtyfold higher, the resonance between the rid — 1 and nj — 2 
dot configurations is now centered at Sd m -4.2 D, and the 
resonance between the = and n</ = 1 dot configurations 
has been reduced to an incipient rise, at the high-ej end of 
the plot. Clearly, the wire potential has displaced the Kondo 
domain towards higher dot-level energies. This displacement 
acknowledged, we recognize in Fig. [4] the salient features of 
Fig. El 

The two landscapes are similar because the phase-shift dif- 
ference on the right-hand side of Eq. (0 is weakly dependent 
on pW. With 5 - 6 W « 7r/2, the conductance curve G(T/Tk) is 
approximately mapped onto G s (T/Tx) throughout the Kondo 
domain. The rise from the high-temperature valley to the 
Kondo plateau is therefore close to universal, dependent on 
the model parameters only through the Kondo temperature 
Tk- The latter is a function of the antiferromagnetic interac- 
tion / between the dot moment and the conduction electrons 
around it. The Schrieffer- Wolff transformation 36 relates that 
interaction to the dot excitation energies: 




Here, Ao (A2) is the energy needed to remove (add) and elec- 
tron to the singly-occupied dot level. For the symmetric 
Hamiltonian (O, in particular, A = A2 = U/2. For nearly 
symmetric Hamiltonians, A = \s^\ and A2 = U + s^. As 
\W\ or \sd + U/2\ grow, the resulting particle-hole asymmetry 



renormalizes the dot energy^^ so that Ao and A2 are changed 
to A* = \s* d \ and A* = U + e*, respectively, where e* is the ef- 
fective dot energy at the LM. 35 

Since both landscapes were computed for the same effective 
width r„, = 0.15 D, only (i) the excitation energies AJ and 
AJ; and (ii) irrelevant operators make the Kondo temperatures 
in Fig. [3] different from those in Fig. [4] The renormalized 
excitation energies displace the Kondo domain along the e ( ; 
axis, while the modified irrelevant operators extend the Kondo 
plateau towards higher temperatures. 

This concludes our overview of the numerically computed 
conductances. Section [VCl will inspect in more detail the data 
in four slices of Figs.[3]and|4]and compare them to Eq. ([9). 



C. Thermal dependence of the conductance 

Figure|5]displays the conductance as a function of the tem- 
perature for U = 5D, r w = 0.15 D, e d + U/2 = 0, and five 
wire potentials: pW = and 1, already studied in Figs.[3]and 
|U and three intermediate values, pW = 0.25, 0.5, and 0.75. 
With W = 0, the open circles represent the symmetric Hamil- 
tonian (0, and the solid line through them reproduces the uni- 
versal function G s (T/Tk) computed in Ref.|2j. Notwithstand- 
ing the wire potentials, the Hamiltonians represented by the 
squares, triangles, and diamonds lie deep inside the Kondo 
regime. For each of them, the phase-shift difference 5 - 6 W in 
Table U is close to n/2. It follows that the right-hand side of 
Eq. (0 is close to G s (T/Tk). The agreement with the numer- 
ical data is excellent. 



7 



1.0 






Ayini| i i iiiiiii 


i i iiiiiii i 


1 IIIIIII 


i iinii| i i nun 

o pW = 0.00 


G 
Q-i 












• pW = 0.25 
□ pW = 0.50 
A pW — 0.75 


0.6 












A rtW — 1 ("If! = 

U = 5D 


0.4 












T„ = 0.15 D 
e d = -2.50 D 


0.2 





















i 


i i 


Mmm"** 




10 


-8 




lO" 6 


10" 4 


T/T K 


lO" 2 



FIG. 5: (Color online) Thermal dependence of the conductance for 
Sa + U/2 = 0, and the indicated values of the other model param- 
eters. The circles, open and filled squares, triangles, and diamonds 
are NRG data, while the solid lines through them depict Eq. ((9}, with 
the Kondo temperatures and phase shifts listed in Table [I] The curve 
through the open circles, in particular, is the universal conductance 
G S (T/T K ) for the symmetric Hamiltonian (O. 2 Since S - <5„. a n/2, 
each solid line is close to G S (T/T K ). 



As the Hamiltonian moves away from the Ed - —U/2 plane, 
the particle-hole asymmetry becomes more pronounced. One 
might expect the difference S - 6 W to grow. As Fig. [2] showed, 
however, in the Kondo regime the Friedel sum rule restrains 
the growth, so that 6 - S w ~ n/2. Illustrative results appear 
in Fig. [6] which displays conductance curves for Sd = -3.4 D. 
Even for the strongest wire potential in the legend, pW = 1, 
the difference 6 - 6 W in Table|I]is only 6% away from n/2. As 
in Fig. [5] therefore, the conductance curves computed from 
Eq. (0 are close to G S {T /Tk)- The agreement with the nu- 
merical data is again flawless. Since we are now closer to 
the boundary of the Kondo regime, the Kondo temperature is 
more sensitive to the renormalization of the dot-level energy 
induced by strong wire potentials. Compared to Fig. [5] Fig. [6] 
thus exhibits a substantially broader spread of crossover tem- 
peratures. 

Figure|7]displays numerical results for Sd + U/2 = - 1 .5 D, a 
still larger departure from the symmetric condition. ForpW < 
0.5, the agreement with Eq. (0 is excellent; for pW = 0.75, 
it is imperfect only at the highest temperatures shown. For 
pW - 1, however, there is substantial disagreement, which 
justifies a digression. 

Inspection of Fig.|4]shows that, for Ed - -AD (and pW - 
1), the model Hamiltonian lies well within the mixed-valence 
regime. 39 In the Kondo regime, Eq. (0 is reliable for thermal 
energies that are small on the scale of the dominant character- 
istic energy Ec — min(|e^|, U + s* d , D). If s, had its bare value, 
Ed = -4 D, the mapping would be reliable for k B T «: E c = D. 
The dot energy has been renormalized, however, and the 
renormalization has pushed the model Hamiltonian into the 
mixed- valence regime. The dominant characteristic energy 
has therefore been changed to E c = min(r w , D) = F w , a reduc- 



Figure Symbol pW 8 w /n S/n 2(6 - 6 w )/n k B T K /D 

5 
5 
5 
5 
5 
6 
6 
6 
6 
6 
7 
7 
7 
7 
7 
8 
8 
8 
8 
8 



TABLE I: Phase shifts and Kondo temperatures for the twenty NRG 
runs depicted in Figs. 15181 The ground-state phase shifts 6 were ob- 
tained from Eq. d!5t , the wire phase shifts 6 W , from Eq. dlOl l. and 
the Kondo temperatures, from the definition G(T = T K ) = Q%/2. As 
explained in the text, a different procedure identified the Kondo tem- 
perature marked with an asterisk, which is associated with a Hamil- 
tonian in the mixed-valence regime. 



tion that restricts the domain of the mapping to kgT <K 0. 15 D. 
The mapping fails at higher temperatures because irrelevant 
operators, which are sizable near the characteristic energy, 
make a significant contribution AG; 1T to the conductance. At 
kgT = 0.1D = 2/3F,,,, for example, the diamonds in Fig.|7]are 
displaced Q.2e 2 /h below the solid line; upon cooling, AG m 
decays in proportion to ksT and becomes insignificant below 
k B T=lQ- 2 D. 

If -Ed were steadily increased beyond -Ed = 4D, the 
model Hamiltonian would traverse the mixed-valence region. 
Once \s* d + U\ > F vv , the dot occupation would approach 
rid = 2. The dominant characteristic energy E c = \s*, + U\ 
would then define the crossover energy scale, which would 
hence rise with -Ed- Soon, the model Hamiltonian would be 
driven to the frozen-level fixed point at the first steps of the 
renormalization-group flow, and the mapping would be re- 
duced to its FL limit, G(T -> 0) = sin 2 (<5 - 5 W ) ~ 0. 

Although Eq. (0 is asymptotically exact at low tempera- 
tures, i. e., for kgT <k E c , as our digression showed its prac- 
tical value is eroded outside the Kondo regime. In the mixed- 
valence regime, in particular, the asymptotic region lies below 
the crossover temperature, i. e., in the vicinity of the FL. To 
plot the rightmost solid line in Fig. [71 we thus had to match 
the right-hand side of Eq. (0 to the diamond at G = 0.702, 
because the identification G(T = Tk) = 0.5^2, which defined 
Tk for all the other plots in Figs. |5][8] became unreliable for 
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FIG. 6: (Color online) Thermal dependence of the conductance for 
Ed + U/2 = -0.9 D. The symbols and lines were computed as de- 
scribed by the caption of Fig. [5] As Table [I] shows, the argument 
6 - S w on the right-hand side of Eq. (|9) is close to n/2. As a con- 
sequence, the solid lines are only slightly different from G s (T/T K ). 
The agreement with the numerical data is, again, excellent. 
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FIG. 7: (Color online) Thermal dependence of the conductance for 
Ed + U/2 = -1.5 D. The symbols and lines were calculated as de- 
scribed by the caption of Fig. [5] As discussed in the text, the high- 
temperature separation between the solid line and the diamonds fin- 
gers a Hamiltonian outside the Kondo regime. 

pW = 1. The asterisk in Table Q] marks the resulting Kondo 
temperature. 

Near the opposite extreme of the Kondo regime, for fixed, 
small -Sd, the wire potential drives the model Hamiltonian to- 
ward the center of the Kondo regime. In Fig. |4] for instance, 
the mixed-valence domain is missing, because the wire poten- 
tial displaced it to positive dot-level energies. Figure [8] dis- 
plays the Sd = -0.4 D plane for the five potentials pW = 0, 
0.25, 0.5, 0.75, and 1. The pW = Hamiltonian is now at 
the boundary of the Kondo regime, and for k B T > 10~ 2 D, ir- 
relevant operators introduce significant deviations AG m from 
the solid line. As pW grows, however, the model Hamiltonian 



sinks deeper into the Kondo regime, and the agreement with 
the solid lines is recovered. 
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FIG. 8: (Color online) Thermal dependence of the conductance for 
Ed + U/2 = 2.1 D. The lines and symbols were computed as described 
by the caption of Fig. [5] The relatively large separation from the sym- 
metric condition E d + U/2 = places the W = Hamiltonian close 
to the border of the Kondo regime; at high temperatures, relatively 
large irrelevant operators, whose influence decays in proportion to 
k B T/D, introduce deviations from Eq. l[9}. Since the wire potential 
displaces the Kondo regime to higher dot-level energies, the distance 
to the border grows with pW. As a result, the squares, triangles, and 
diamonds are in excellent agreement with Eq. (|9). 



VI. CONCLUSIONS 

In the Kondo regime, the physical properties for the Ander- 
son model are universal functions of the temperature scaled 
by the Kondo temperature. The conductance constitutes no 
exception. For e,/ + U/2 — W = 0, the Anderson Hamiltonian 
reduces to Eq. (0, in which case the thermal dependence of 
the SET conductance is the universal function G s (T /Tk)^ 
For nonzero Ed + U/2 or W, Eq. (0 maps linearly the conduc- 
tance onto G S (T/T K )M- With e d + U/2 = W = 0, the particle- 
hole symmetry of the Hamiltonian © forces the ground-state 
phase shift 5 to be n/2 and reduces Eq. (0 to a trivial identity. 
For asymmetric Hamiltonians in the Kondo regime, as the 5th 
and 6th columns in Table Q] and the plots in Fig. [2] suggest, 
while the ground-state phase shift can take any value in its do- 
main of definition, the Friedel sum rule keeps the difference 
6 - 6 W close to n/2. The linear coefficient of the mapping (0 
is never far from -1, and the thermal dependence of the SET 
conductance, never far from G s (T/Tk). 

Equation (0 becomes asymptotically exact as k B T <sz E c . 
In the Kondo regime, the dominant characteristic energy is 
E c = min(Z), |e^|, U + s* d ), and the mapping to the conductance 
is reliable throughout the crossover from the LM to the FL. In 
the mixed-valence regime, with |e^| < r u , or U + s*. < T w , the 
dominant characteristic energy E c = min(Z), F H ) reduces the 
domain of the universal mapping to a temperature range close 
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to the low-temperature fixed point, i. e., to the final steps in 
the rise to the limit G(T = 0) = Q 2 sin 2 (5 - 5 W )£- 



Summary 

To recapitulate, the essentially exact numerical data in 
Figs. I23IH] offer an overview of electrical conduction through 
a quantum dot embedded in the conductance path of a nanos- 
tructured device. The linear mapping (O to the universal func- 
tion G s (T/Tk) for the symmetric Anderson Hamiltonian 2 ' 8 
describes accurately the thermal dependence of the conduc- 
tance for IcbT <sc D in the Kondo regime. A gate potential ap- 
plied to the wires affects only quantitatively the dependence 
of the conductance on the temperature and dot-level energy. 



In the Kondo regime, independently of the wire poten- 
tial, the Friedel sum rule drives the thermal dependence of 
the conductance to the neighborhood of G s (T/Tk)- In par- 
ticular, G(T <K Tk) ~ 02- In the mixed-valence regime, 
only well below the crossover temperature does G(T/Tk) map 
onto G s (T/Tk), and G(T — > 0) is substantially smaller than 
the quantum conductance. In the laboratory, it is often diffi- 
cult to distinguish the Kondo regime from the mixed-valence 
regime. 19,20 The mapping to the universal curve offers a prac- 
tical solution to this problem. 
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